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defined by a cylindrical detector and a point 
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' We derive analytical expressions for the solid angle subtended by a right circular 

■ cylinder at a point source with cosine angular distribution in the case where the 

, source and the cylinder axes are mutually orthogonal. 

o 
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The calculation of the solid angle subtended by a cylindrical detector at a 
^ point source is of common interest in nuclear science. The case of an isotropic 

source has been extensively studied, but, to the best of our knowledge, no 
work has been published considering a point cosine source. Such source could 
arise as the second term in a Legendre expansion of a general source or, to take 
a well known example, as one of the terms in the Fermi expression = 
1 + \^3fx. This expression holds for the angular distribution of low energy 
neutrons leaking from a variety of scattering materials and, if suitabl y modified 



1 Introduction 



= 1 + Aji) , can also be used in the case of lithium hydride (jVerbinski 



1967l b The scattering of neutrons from hydrogen nuclei considered at rest is 



also described by a cosine distribution 



In this work the situation of a point cosine source defined with respect to some 
axis and a right circular cylinder with axis orthogonal to that of the source 



Email address: mjprata@sapo.pt (M. J. Prata). 
1 Supported by Fundagao para a Ciencia e Tecnologia (Programa Praxis XXI 
BD/15808/98) 



Preprint submitted to Elsevier Science 



7 February 2008 



is considered. Under this restriction we present the analytical calculation of 
the solid angle subtended at the source positioned at an arbitrary location. 
Sample graphics of the expressions obtained are presented. 



2 Solid Angle Calculation 



Let the unit vector k define the source direction. The source distribution 
giving the probability f(fl)dfl of emission around the direction of the 
unit vector ft is defined by 



/(«) 



^ ,k- n > o 

IT ' 

o ,k- n < o 



The distribution is normalized so that 



jj f(n)<m = i. (2) 

all directions 

The solid angle (Q) is given by 

n= II f(fi)dn. (3) 

directions 
hitting detector 

The origin of the coordinate system is chosen to coincide with the point source, 
the x axis aligned with k and the z axis parallel to the cylinder axis. In figs. 1 to 
3 we show the three cases to be considered and introduce some of the notation 
used. Generally the solid angle is a sum of two terms, the one subtended by 
the cylindrical surface (fl cy i) and the other by either of the end circles (Q circ ). 
In case (i) one has Q = Vt cyh in case (ii) Q = Q cy[ + Vt circ and in case (iii), 
n = Q C i rc . In the various cases L 1; L 2 , d and r are all positive. From the 
symmetry of the problem, the solid angle is an even function of a and in the 
following we will thus consider a > 0. In cases (i) and (ii) one has d > r and, 
in case (iii), r > d. With the notation used there results that 

k ft = cos(a + ip) sin(0) (4) 

where 8 is the polar angle from the z axis and tp is the azimuthal angle in the 
xy plane as measured from an axis through the origin and intersecting the 



2 



detector axis. Hence 




Fig. 1. Geometry in Case (i) 




Fig. 2. Geometry in case (ii) 



2.1 Integration limits 

Because the source only emits into the hemisphere corresponding to i > 0, 
the range of variation of ip depends on the value of a, as explained in figs. 4 
to 7, where the the hashed area shows the illuminated part of the detector. 
Introducing 

ipQ = arcsin (r/d) , (6) 

there results 

ai = 7r/2-y?o (7) 
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and 



Fig. 3. Geometry in case (iii) 
a c = 7r/2 + ip . 



(8) 
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Fig. 4. Definitions of ipo, a\ and a c (cases (i) and (ii) ) 
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Fig. 5. Integration limits for ip in cases (i) and (ii) when < a < ot\ 
The limits for tp are then given in table 1. 

The limits for 9 can be expressed in terms of p_((p) and p+(p) that were 
introduced in the pictures and are given by 

p ± ((p) = dcos(p) ± \Jr 2 — (dsm(p)) 2 . (9) 
4 



Fig. 6. Integration limits for ip in cases (i) and (ii) when ot\ < a < ot c 




Fig. 7. Integration limits for ip in case (iii) 



Table 1 

Integration limits for <p 



Case a range ip min ip max 

< q < cci —tpo ipo 
(i) or (ii) «i < q < a c —ipo ir/2 — a 
a c < a < T[ = 

(iii) < a < it —it/2 — a it/2 — a 




Fig. 8. Angles 4>± in cases (i) and (ii) (d ^ r ) 
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Fig. 9. Angle 4> + in case (iii) (d < r) 



Table 2 

Integration limits for 9 



Case £l cy i or Q, 



(i) 




(ii) 


^cyl 


(ii) 




(iii) 


Qcirc 



vr/2 - arctan(^) 
vr/2 + arctan(|i) 
vr/2 + arctan(^) 
vr/2 + arctan(^) 



vr/2 + arctan(^-) 
vr/2 + arctan(|2.) 
vr/2 + arctan(^) 

7T 



The limits for 6 are summarized in table 2. 



2.2 Integration 



Let us define 



fi(L,d,r,a,ip) = /i(</?) 



2tt 



cos(a + </?) 



arctan 




L 2 + p 2 



dtp 
(10) 



and 



f 2 (L, d, r, a, ip) = f 2 ((p) = J cos(a + ip) 



arctan — + 



L \ Lp + 



P+ V + pl 



dip 
(11) 



Then, performing the first integration in the rhs of eq. 5, where # min and E 
are given in table 2, yields in each case: 
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case (i) 

n = n cy i = lf 1 (L 1 ,d,r,a, V ) + f 1 (L 2 ,d,r,a, V )f™: , (12) 

case (ii) 

tocyi = [fi(L 2 , d, r, a, <p) - d, r, a, J™ , (13) 

n circ =\f 1 (L 1 ,d,r,a,<p)- f^d^a,?)]^ , (14) 

fi = fi cy , + fi cire = [/i(L2,d,r,a,^)-/2(Li,d,r,a,^)]J^ , (15) 

case (iii) 

Q = Q circ = - cos(a + v?)^-/ 2 (^i,rf,r,a,^)|^ n x , (16) 

= i-/ 2 (L 1 ,d,r,a, ¥ »)fc. (17) 

Eq. 17 is obtained from eq. 16 by integration of the 1st term using the values 
from table 1. 

Using again the limits from table 1, eqs. 12 to 15 and 17 can be rewritten as 
case (i) 

n = Vl cyl = Fi(Li, d, r, a) + Fx(L 2 , d, r, a) , (18) 

case (ii) 

Qcyi = F 1 (L 2 ,d,r,a) - Fx(Lx, d, r, a) , (19) 
^circ = Fi(Li,rf,r,a) -F 2 (L 1 ,d,r, a) , (20) 

Q = Q cyi + Q circ = F 1 (L 2 , d, r, a) - F 2 (L ± , d, r, a) , (21) 
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Table 3 

Definition of the integrals 



Integral a range Value ¥> m in ¥>n 



< a < ai ML, d, r, a, ip)\^Z -ipo ip 



(i = l,2) a!<a<a c /;(L, d, r, a, <p)\^ -<po n/2-a 



Ct c < Ct < IT 



F 3 0<a<7r 1/2- f 2 (L,d,r,a,p)\^ -vr/2 - a vr/2 - a 

case (iii) 

Q = n circ = F 3 (L 1 ,d,r,a) , (22) 
where the integrals Fj (i = 1..3) are defined in table 3. 

,2. ,2.1 Calculation of Fx 

Integrating by parts the 1st term in the rhs of eq. 10, 

L \ , , , ( L 



J cos(o; + ip) arctan y — J dip = sin(o; + ip) arctan y — J 

+ / sin(a + y , ™~ Jbp (23) 



Substituting eq. 23 in the rhs of eq. 10 yields, after some algebra, 

f x {ip) = — [A w + cos(a)A u + sm(a)A 12 \ , (24) 

Z7T 



where 



Aio = sin(a + </?) arctan f — ^ , (25) 



^-/toH^t^&f 1 * (26) 
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and 



/ 



Lp_ I p_ sm(ip) 



L 2 + P 2 _ \^r 2 - [dsin(p)] 2 ^ ' 
Performing a change of integration variable to 0_ represented in fig. 8, 



(27) 



2 arctan 



P- sin(^) 



d + r — p cos(v?) J 



(28) 



we find 



COs(0_ 



i? + pL(4>-) 



d4>. 



(29) 



and 



A 12 = Lr 



sin 



L 2 + p 2 _(0_) 



(30) 



where 



p_((f)-) = \J 1 d 2 + r 2 — 2dr cos(0_) 
One then easily obtains 



(31) 



* 1 = 53 



a/1 — m 2 



arctan | \l — — — tan ( — 

1 — m V 2 , 



(32) 



where 



2dr 



L 2 + d 2 + r 2 



; < m < 1 



(33) 



and 



A 2 = 2rflog(L 2 + p 2 ). 



(34) 
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2. 2. 2 Evaluation of F 1 (d^r) 

For < a < a± the terms proportional to sin (a) vanish and one gets 



+ 



d \ a/1 - ^ 



Fi(L,ci,r,a) 



arctan 



cos(a) 



7T 



d 



arctan 



Vd 



1 + m I d — r \ I Id — r 

— arctan 



1 — mV d + r 



d + r 



(35) 



For «i < a < a c , the expressions are less simple and we give each term in eq. 
24 separately. Setting 



P- = P-(y?max) = P_\<p=n/2-a 



= dsin (a) — yr 2 — [dcos (a)] , 
the first term can be written as 



(36) 
(37) 



^4io|- max — arctan | — | — arctan 



' V^min 



L 



sin(a - (p ) 



Defining 



and 



a\ = tan 



'd — r 
d + r 



(38) 
(39) 



bi = tan 



cos (a) p_ 



then 



V^max 



r + d — sin (a) p_ 



(40) 



A I (Pmax 



L 

2d 1 VI -m 2 



arctan 



' 1 + m 
1 — m 



&i — arctan 



-2 [arctan (frx) — arctan (cii)]} 



Finally, from results 



1 + m 
1 — m 



ai 



(41) 



A 12 |£- = illog[- 



2d OL L 2 + rf 2 -r 2J 



(42) 
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Using eqs. 24, 38, 41 and 42, F\ can be calculated in the case a\ < a < a c , 

2. 2. 3 Evaluation of F 1 (d = r) 

The value of F 1 when d = r can be obtained from eq. 10 using the integration 
limits from table 3 or calculating the limit of the expressions obtained in 2.2.2, 
for a > cti. When L^O one obtains Fi(L, r, r, a) = (l + cos«)/4. Since Fi = 
for L = 0, Fi is discontinuous when r = d, L — > 0. 

2.2.4 Calculation of F 2 

The calculation is very similar to that of F±. Integration by parts gives 




dip = sin (a + ip) arctan 





(43) 



Eq. 11 can then be rewritten as 



f2((p) = — [A 20 + cos(a)A 21 + sm(a)A 22 \ , 



(44) 



where 




(45) 




(46) 



and 




(47) 



The integration variable is now changed to <p + (fig. 8) 




(48) 
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and the integrals are then expressed as 



A 2 i = Lr 



I 



cosQ+) 
L 2 + p\ (</>+) 



(49) 



A 22 = -Lr 



I 



sin 



L 2 + p\ (0 + ) 



d(f>+ , 



where 



(50) 



p + (0+) = \Jd 2 + r 2 + 2dr cos(0_| 
The integration is straightforward resulting in 



(51) 



A * = 2d 



2 / ll-m 

arctan \ / tan — 

l + m \ 2 



a/1 — m? 



and 



^22 = ^ l °^ L2 + -°+) • 



(52) 



(53) 



where m is obtained from eq. 33. 
2.2.5 Evaluation of F 2 (d ^ r) 

For < a < a± the terms proportional to sin (a) again vanish giving 



d \ a/T 



F 2 {L,d,r,a) 



: arctan 



cos(a) 



7T 



r L 
— arctan = 
d \^d 2 -r 2 



1 — m d + r 
1+mVd-r 



arctan 



>d + r 
d— r 



(54) 



It is interesting that using arctan(z) + arctan(l/z) = n/2 the expression can 
be recast as 



F 2 (L, d, r, a) = F^L, d, r, a) + cos(a) L/(2d) (l - - m 2 ) . (55) 

For a\ < a < a c , setting 
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P+ — P+( ( / 9 max) — P+\ifi=Tr/2-a 



= dsm (a) + \Jr 2 — (dcos (a)) , 
the first term is determined from 



^2olw max = arctan f — ) — arctan 



Vd 2 ^ 



Considering 



a 2 = tan 
b 2 = tan 



<d + r 
d — r ' 
cos (a) Ph 



'/'max 



r — d + sin (a) p + 



gives 



(56) 
(57) 



sin(a — <f ) . (58) 

(59) 
(60) 



-^21 1- 



arctan 



2d 



{2 [arctan (6 2 ) — arctan (a 2 )] 



1-m. | I 1 — m 

— 2 - arctan W — a 2 

1 + m / i V 1 + m 



(61) 



Finally, for A 22 the result is 



A 22 \^ 



2d ° g[ L 2 + d 2 -r 21 



which completes the evaluation of F 2 when r ^ d. 



(62) 



2.2.6 Evaluation of F 2 (d — r) 

Taking the limit of the expressions obtained in 2.2.5 when L ^ 0, one obtains 
F 2 (L,r,r,a = 0) = 1/2 + L/(2r)(l - - m 2 ), F 2 (L,r,r,a = tt) = and 



F 2 (L, r, r, a) = (1 + cos a)/4 — 1/(2tt) {arctan [2(r/L) sin a] 
+ cos(a)L/r (a — 7r + 1/Vl — m 2 (ix /2 + arctan cot(a)Vl — m/y/l + m 

- sin(a)L/(2r) log [l + 4r 2 sin 2 (a)/L 2 l } . (63) 
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for sin(a) ^ 0. A simple calculation shows that when L — > 0, F 2 (L, r, r, a) — > 
0, regardless of a This is also true when d ^ r, and F 2 is thus continuous 
when L — > 0. 



2.2.7 Calculation of F 3 
From table 3, 



F 3 (L, d, r, a) = 1/2 - / 2 (L, d, r, a, ^)|^ . (64) 

The calculation of / 2 has already been done although the quantities involved 
(0 + , p + ) have a different graphical interpretation as shown in fig. 9. Eqs. 48 
and 52 can still be used but some care is required as explained further on. 
From eqs. 64 and 44 follows that 

Fs = h [Aso + cos ^ Asi + sin ^ A ^z: > ( 65 ) 

where 



A 30 = ip - A 20 , (66) 
A 31 = -A 21 , (67) 
A 32 = -A 22 (68) 



and A 2i , (i = 1..3) are determined from eqs. 45, 52 and 53. 

2. 2. 8 Evaluation of F 3 (d ^ 0) 

With y? max and </? m i n chosen from table 1 



P+(y?max) = P+ , (69) 

P+(Win) = -P- , (70) 

where p_ and p + are found from eqs. 37 and 57. To calculate ^3i|SX the 
quantities 



a 3 = tan 



cos (a) p_ 
r — d + sin (a) p_ 



(71) 
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and 



b-i = tan 



will be useful. A 30 is given by 



V-'max 



cos (a) p + 
r — d + sin (a) p + 



(72) 



A 3o|^ n x = 7T - [arctan (L/p + ) - arctan 
Using eq. A.l from the appendix results in 



(73) 



4 3 n|^ max = vr - 2 arctan (l/z) , 



where 



2; = 



r 2 - d 2 - L 2 + yV 2 - d 2 - L 2 ) 2 + 4L 2 (r 2 - [dcos(a)] 2 ) 



2L^/r 2 - [rfcos(a)] 2 
Then, applying arctan (z) + arctan (l/z) = tt/2 yields 



-430 1 - max = 2 arctan (z) . 



(74) 



(75) 



(76) 



Some care is required to obtain A 3 \. From fig. 9 is clear that as a — > 7r/2 then 
0+ (Vmin) - ► — 7T and that — 7r < + (</? m i n ) < Sir/ 2 for ct > it/ 2. Therefore 
a 3 = tan [0 + ((/9 min )/2] is not continuous as a goes through tt/2. The discon- 
tinuity also appears in the rhs of eq. 48 when ip decreases to less than —it. 
From eqs. 67 and 52 its obvious that both terms in A 31 will be discontinuous. 
The integral, of course, must be continuous and a proper expression can be 
provided by noticing that while 



z 7^ arctan [tan(z)] ; — it < z < —tt/2 , 



one has 



(77) 



z = —tt/2 — arctan [1/ tan [z)\ ; — it < z < 
Thus, to evaluate A 31 , the substitutions 



(78) 



arctan(a 3 ) — > — tt/2 — arctan(l/a 3 ) , 



(79) 



arctan 



'1 — m 
1 + m 



a 3 



-tt/2 — arctan 



(80) 
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are used and it is found that 



A I '/'max 

A 3l| Vmin 



(81) 



where 



At i = arctan(6 3 ) + arctan(l/a 3 ) , 



(82) 



M — j yi M _|_ j Yt 

At 2 = arctan(W b 3 ) + arctan(W /a 3 ) 

Vl + m V 1 — m 



(83) 



The sum of arctangents can be cast as a single one using eq. A.l from the 
appendix and after some algebra one gets 



At i = 2 arctan [d cos (a) / (r + E)] , 



(84) 



At 2 = 2 arctan 



{Gd cos (a)} I \EE + r^J H 2 - [2Ldcos (a) 



(85) 



where 



E = yr 2 — [dcos(a)] , 



(86) 



G = L 2 + d 2 - r 2 



and 



H 



L 2 + (d + r) 2 } \L 2 + {d-rf 



The remaining term is given by 



(87) 



A I ^max — 



2d g [L 2 + fP 



(89) 



5. 2. 9 Evaluation of F 3 (d = 0) 



Calculation of the limit of the expressions obtained in 2.2.8 yields F 3 (L,d = 
0, r, a) — 1/2 — (arctan [L/r] + Lr/ (L 2 + r 2 )). In the same way one can show 
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that F 3 (L — > 0, d < r,r,a) = 1/2 regardless of d, r and a. Altough F 3 is 
only needed fo 

1/4(1 + cos a). 



only needed for d < r, it is worth mentioning that lim( lim F 3 (L, d, r, a)) = 

L— >0 d- ~ 



3 Sample graphics and discussion 



To illustrate the behavior with respect to a we consider a cylinder of lenght 
10 and radius 1. In figs. 10, 11 and 12 are shown examples of cases (i),(ii) and 
(iii), respectively. As argued before, Q is an even function of a in all cases. In 
cases (i) and (ii), for0<|a|<7r/2 — arcsin(r/<f), the solid angle is simply 
proportional to cos(a); for \a\ > n/2 + arcsin(r/rf), Q = 0. In between, the 
dependence on a is more complicated because a fraction of the cylinder is 
not illuminated by the source, but Q(a) obviously decreases as |a| increases. 
The transition region is large for d m r (d < r) and becomes increasingly 
narrower as d increases. For d ^> r, the region essentially vanishes so that 
Q oc cos(a) when \a\ ;$ n/2 and £1 = when \a\ <J vr/2. It is intuitive that 
in case (i), Q decreases when distance d increases and all other parameters 
are held constant. This is not necessarily true for case (ii), because source 
and cylinder are in a skew geometry. Considering two distances d\ < d 2 , one 
concludes from the preceding discussion that there is certainly some region for 
a where Q(d 2 ) = and Q(di) > so that Q(di) > 0(cZ 2 ) ■ As |a| decreases to 0, 
this relation can be inverted, which happens for all distances shown in fig. 11. 
Since for very large d, £1 is certainly a decreasing function of d, there must be 
a maximum, for each a. This is illustrated in fig. 13, for a = 0, L x = 5 and two 
cylinder lengths: 10 (L 2 = 15) and 20 (L 2 = 25). It is worth mentioning that 
in case (ii), apart form the a dependence, a similar behavior can be observed 
for an isotropic point source (i.e. a maximum when d is changed and L\ and 
L 2 kept constant). Case (iii) is exemplified in fig. 12. In this case (r > d) the 
solid angle is defined only by one the end circles which is, in general, partially 
illuminated by the source. The exception happens when d — ► r : for a = 
the circle is fully illuminated where as for a = ir it is totally obscured. When 
d = the source is aligned with the center of circle and Q is then independent 
of a. As d increases, for L\ ^ 0, the dependence on a becomes stronger with 
a maximum when a — 0. When L\ — (not shown in fig. 12), Q = 1/2, 
regardless of a or d. 
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Fig. 10. Solid angle in case (i), for a cylinder with radius 1 and length 10 = 5 + 5 
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Fig. 11. Solid angle in case (ii) for a cylinder with radius 1 and length 10 = 15 — 5 
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Fig. 12. Solid angle in case (hi) when O = £l c irc, for a circle with radius 1. The 
intersection of the three curves near ir/2 is not exact. 



4 Conclusions and outlook 



The solid angle denned by a point cosine source and a right circular cylinder 
with axis orthogonal to that of the source has been treated analitically. It 
has been shown that for d > r and ^ |a| < ir/2 — arcsin (r/<f), the whole 
dependence on a is given by a cos(a) factor (eqs. 35 and 54). It is possible 
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0.025 
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0.01 
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° 5 10 15 20 

d 

Fig. 13. Peak values (a = 0) of the solid angle in case (ii) as a function of d, for two 
cylinders of radius 1 and lengths 10 (15 — 5) and 20 (25 — 5) 

to obtain, in a similar fashion, analytical expressions in the case of a cosine 
source distributed on a wire paralell to the cylinder axis. A work where we 
report these results is in preparation. 
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A Sum of Arctans 



To prove the identity 

arctan x ± arctan y = 2 arctan (x ± y) / ^1 =F X U + \j {x ± y) 2 + (1 =F xy) 2 ^j 

(A.l 

consider h(x, y) = arctan a; ± arctan?/ and g(x, y) = 2 arctan [(x ±y) / D] where 



D = (l t xy + ^/{x±yf + {lTxyf^j 



Using tan(w±f ) = [tan(-u) ± tan(v)] / [1 =F tan(w) tan(w)], we find that tan (h) = 
(x ± y) / (1 =F X U)- I n a similar way, applying tan(2w) = 2 tan(-u)/ [1 — tan 2 (w)], 
gives 



tan [g{x,y)\ = 2 



(x±y)/D 



D 2 - (x± y y\ /D 2 
(x±y)/(l^xy) , 
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where in the last step we used the fact that D ^ 0. Since tan (h) = tan (g) 
then h — g = nir for some n. Taking x = y = results that n = and, because 
h and g are continuous, we conclude that eq. A.l holds for any x and y. 
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